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Abstract 

We propose a two band model for superconductivity. It turns out that the 
simplest nontrivial case considers solely interband scattering, and both bands 
can be modeled as symmetric (around the Fermi level) and flat, thus each 
band is completely characterized by its half-band width Wn (n=l,2). A use- 
ful dimensionless parameter is 6, proportional to W2 — Wi. The case 5 = 
retrieves the conventional BCS model. We probe the specific heat, the ratio 
gap over critical temperature, the thermodynamic critical field and tunneling 
conductance as functions of 6 and temperature (from zero to Tc). We com- 
pare our results with experimental results for MgB2 and good quantitative 
agreement is obtained, indicating the relevance of interband coupling. Work 
in progress also considers the inclusion of band hybridization and general 
interband as well as intra-band scattering mechanisms. 
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I. INTRODUCTION 



Magnesium Diboride [MgB2) appears to be a rather "unconventional" conventional 
superconductor^'^. Two band effects observed as deviations of conventional BCS include: 
anomalous specific heat^ and two gaps features (including double peaked tunneling conduc- 
tance spectra)^^^*^. The superconductive mechanism, nevertheless seems to be conventional 
phonon BCS-likc^^. In this short communication we present a two band model based on 
the classical work by Suhl et al.^^ and on an extension of the latter applied to high Tc 
compounds^^. We mention other multiband models in the literature^^"^^, and some calcu- 
lations and fittings within a multiband and strong coupling context include Ref. 22-25. In 
section II we introduce a two band modeP^'^^ and within the usual BCS scheme we compute 
the mean field expressions for the free energy, entropy, critical field, conductance, and the 
selfconsistent equations for the gaps functions. In particular we consider the simplest case: 
solely interband pairing coupling via phonons. In section III we compare our simple model 
with some experimental results for the case of MgB2^^'^''^'^ , indicating that the interband 
pairing mechanism is somehow relevant. Finally in section IV we present some concluding 
remarks and future work. 



II. THE TWO BAND MODEL 

Our model follows Ref. 12, 13, with the Hamiltonian 

k,m kq,m 

where the c[,'s are the usual creation operators, E^^m ai'e the bands dispersion (m = 1,2), 
Vn^m are the positive pairing coefficients (V12 = V21 and D = V11V22 — V^il 7^ 0). We have 
defined k = (k, |), —k = (— k, j), iV is the number of sites and the last summation is with 
the usual energy cutoff cud ■ The order parameters A„ are defined as the expectation values 

k,m 

The effective Hamiltonian is given by (within the Hartree Fock scheme for anomalous 
pairing , see Ref. 13) 

k,m 

where 

and (Ta;, az are the usual Pauli matrices. The free energy per site F is given by 

exp {-I3NF) = Trexp {-pH^^ 
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k,m 



where f{u) = {exp{(3u) + ujk,m = ^E^m + and fk,m = f{i^k,m)- The relative free 
energy 6F — F — F(Ai = A2 = 0), the thermodynamic critical field entropy (per site) 
and specific heat are given, respectively by 

(1 + cosh PUk,m) _ 1 

k,m 



SF{T) - ^0 - ^ E In (i + eosh/3£;.,J - - 8^^^^ 
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-5 = E ((1 ~ /fc-'") ^n(l - /fc,m) + /fc,mln/fc,m) (3) 

fe.m 



/ 1 dA^ \ 

j:hm(^-h,m)U^ + ^f3^] (4) 
k,m \ / 



The condensation energy is given by 

5F{T ^0) = Wc = Eo-^Y. (^fc-- - ^fc--) (5) 

and the superconductor- normal tunneling differential conductance (conveniently scaled) is 
defined by 

GiV) = - E / depmAe f^'^^^ (6) 



Pm,s{£) = Pm sign{e)Je'^ - Real 



r ^ 0^ 



Minimization of the free energy with respect to the gaps functions, yields a coupled 
nonlinear system of integral equations for the gaps, to be solved self consistently, and given 

by 

(V22 - DR,{Ai, T)) Ai - VuA2 = 

(7) 

-I/12A1 + {Vu - DR2{A2, T)) A2 = 



where 



i?„(A^, T) = dePm{e)S ^e^ + Aij , S{x) = ^ tanh ( 



X 

2r 



and with Pmi^) the density of states associated to the respective band. The transition 
temperature is the highest temperature Tc — solution of 



{V22 - DR,{0,T,)) {Vn - DR2{0,T,)) = V, 
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III. RESULTS 



We compute the observables presented in the previous section. In particular we consider 
only interband scattering Vu = V22 = 0, V12 = A, the simplest relevant case-*^^'-*^^. We consider 
two flat symmetric bands, with Pmi^) = Pm(0) = Pm- 

The gaps equations (7) now read 

= Ap„A„i?(A„, T), n 7^ m = 1, 2 




At zero temperature the gaps equations are given by (in convenient units) 

2Ai(r ^ 0) 1 

2Ai(T = 0) 1 / 1 

02 = — 5^^77^ — = exp - 1 - - 



(8) 



where = ^piP2 and a satisfies a selfconsistent equation. An excellent approximate 
solution is given by 



with 

-1<<5=^^<1 

Pi +P2 

Notice that all the above mentioned observables will yield the standard BCS expressions^^ 
in the limit 5 = 0. 

The critical temperature is given by Tc — l.lSa;^ exp(— ^■~^). 

We label the bands such that 5 > 0. If we consider MgB2, from Ref. 1, 27 we have 
Tc ~ 40 ^K, cun ~ 800°ir yielding ^ ~ 0.32. From Ref. 26 we approximate Wi ^ 5.6eV 
(pi ^ 0.179e1/-i), W2 ^ U.eV (p2 ~ 0.071e\/-i) yielding 6 ^ 0.432. 

In Figure 1 we plot the normalized gaps 0^ at zero temperature, Eq.(8), and (minus) 
the condensation energy —Wc Eq.(5), both as function of S. The condensation energy is 
normalized to the BCS reference state i.e. Wc = SF{S,T = 0)/SF{S = 0,T = 0) (see 
Eq.2). The chosen normalization yields the standard BCS (weak coupling) value of unity 
for the gaps and the condensation energy. As 6 is varied away from zero the condensation 
energy is less than the standard BCS. One gap will depart from weak to 'a medium couphng 
regime' (02 > 1), conversely the other gap will dive towards 'a less than weak coupling 
regime' (0i < 1), with the geometrical average (f)i{S)(f)2{S) = 1 always in the standard weak 
coupling regime. These features seem consistent as we fit the parameter 8 with experimental 
data^-^°. 
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In order to solve for the gaps, Eq.(7 ), we can use the available low temperature and 
near the critical temperature expansions'^. These allow us to nicely interpolate, for the full 
temperature regime < t = T/Tc < 1. Once this is done we can readily compute the 
specific heat, Eq.(4), entropy, Eq.(3), and the thermodynamic critical field, Eq.(2). 

In Figure 2 we plot the specific heat Cy (normalized to the normal state value at Tc) 
versus the temperature r for several values of S. The standard BCS result is represented 
by the curve S = 0. The anomalous behavior of Cy consists in going under the BCS value 
in the region 0.5 < r < 1, and going over the BCS value in the region < r < 0.5. 
This feature is in very good agreement with Ref. 3. In Figure 3 we plot the entropy S 
(normalized to the normal state value at Tc) versus the temperature r for several values of 5. 
The standard BCS result is represented by the curve 5 = 0. As 5 departs from zero (bands 
are less 'identical') the system increases its entropy. In Figure 4 we plot the thermodynamic 
critical field (normalized to the reference state 6 = 0, T = 0) versus the temperature r for 
several values of S. The standard BCS result is again represented by the curve 5 = 0. As 5 
increases the critical field is reduced when compared to the BCS value. This is in agreement 
with experimental results^. 

In Figure 5 we plot the conductance, Eq.(6) versus applied voltage, for a fixed value of 
5 — 0.5, and for several temperatures, and where a small dispersion is included, F = 0.1 
meV'^.The double peaked form is in very good agreement with observations (see for example 
Ref. 7). 



IV. CONCLUDING REMARKS 

We presented the simplest relevant two band model for superconductivity, based on a 
standard BCS-like pairing mechanism. We computed the gaps equations at zero tempera- 
ture. Also the specific heat, entropy, critical field and conductance as function of tempera- 
ture. We considered the simplest interband scattering mechanism (one pairing parameter) 
and two planar symmetrical bands (one parameter band model). Our results seems to 
be in very good agreement with some experimental results on the compound MgB2, in- 
dicating that interband pairing is somehow relevant for this compound. These results are 
being investigated further. Work in progress incorporates intraband pairing mechanisms, an 
hybrid- like interband pairing mechanism^^, absent in most of the theoretical models, and a 
more involved band structure. 
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Figure Captions 



Figure 1: Gaps 0^ at zero temperature and (minus) the condensation energy —Wc versus 
S. Convenient units 0^ = 2A^/3.537; and = SF{S, T = 0)/SF{S = 0, T = 0). See text. 

Figure 2: Specific heat Cy versus r = T/Tc for several 6 values, normalized to the normal 
state specific heat at T^. CniTc) = 47r^(pi + p2)/6. See text. 

Figure 3: Entropy S versus t = T/Tc for several 5 values; normalized to the normal state 
entropy at Tc. Sn{Tc) = in^ipi + p2)Tc/6. See text. 

Figure 4: Thermodynamic critical field versus r = T/T^ for several 6 values, normalized 
to the reference (BCS) state T = 0, 5 = 0, H^{t) = F{t, 5)/F{0, 0). See text. 

Figures 5: Tunneling conductance for several temperatures, at 5 — 0.5, versus applied 
voltage. A small dispersion is included F = 0.1 meV^^. See text. 
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